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ABSTRACT 

We utilize Sloan Digital Sky Survey imaging and spectroscopy of ^53,000 star-forming galaxies at z ~ 0.1 
to study the relation between stellar mass and gas-phase metallicity. We derive gas-phase oxygen abundances 
and stellar masses using new techniques which make use of the latest stellar evolutionary synthesis and pho- 
toionization models. We find a tight (±0.1 dex) correlation between stellar mass and metallicity spanning 
over 3 orders of magnitude in stellar mass and a factor of 10 in metallicity. The relation is relatively steep 
from 10 8 5 - 10 10 5 Mq/z^q, in good accord with known trends between luminosity and metallicity, but flat- 
tens above 10 10 5 M©. We use indirect estimates of the gas mass based on the Ha luminosity to compare our 
data to predictions from simple closed box chemical evolution models. We show that metal loss is strongly 
anti-correlated with baryonic mass, with low mass dwarf galaxies being 5 times more metal-depleted than L* 
galaxies at z ~ 0. 1 . Evidence for metal depletion is not confined to dwarf galaxies, but is found in galaxies with 
masses as high as 10 10 M . We interpret this as strong evidence both of the ubiquity of galactic winds and of 
their effectiveness in removing metals from galaxy potential wells. 

Subject headings: galaxies:abundances — galaxies: evolution — galaxies: fundamental parameters — galax- 
ies: statistics 



1. INTRODUCTION: 

Stellar mass and metallicity are two of the most fundamen- 
tal physical properties of galaxies. Both are metrics of the 
galaxy evolution process, the former reflecting the amount 
of gas locked up into stars, and the latter reflecting the gas 
reprocessed by stars and any exchange of gas between the 
galaxy and its environment. Understanding how these quanti- 
ties evolve with time and in relation to one another is central 
to understanding the physical processes that govern the effi- 
ciency and timing of star formation in galaxies. 

The influence of stellar winds and supernovae on the 
interstellar medium (ISM) of galaxies, generally dubbed 
'feedback', has been regarded as an imp ortant in gredi- 
ent in galaxy evolution since the 1970s JT arsorJ fl97l 
Larson & Dinersteidll975t IWhite & Reeslll978l) . Feedback 
is believed to play a critical role in regulating star formation 
by reheating the cold ISM and by physically removing gas 
from the disk and possibly the halo via galactic winds. 
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Various lines of observational evidence have established that 
large-scale outflows of gas are ubiquitous among the most 
actively star-forming galaxie s at low and high redshift (e.g. 
Lehnert & Heckman 1 19961 iDah lem. Weaver, & Heckman 
1998HRunke. Veilleux. & Sandersll2002HShaplev et alJl2001t 
iFrve. B roadhurst. & Bemtez 2002). However, detailed 
studies of winds in nearby starbursts have shown them to be 
a complex, multiphase, hydrodynamical phenomenon with 
the majority of the energy and newly synthesized metals 
existing in the har d-to-observe coronal a nd hot phases (T = 
10 5 - 10 7 K) (e.g. Stric kland et alJl2002l) . This complexity 
has prevented both a direct assessment of the cosmological 
impact of galactic winds, and the development of physi- 
cally accurate prescriptions for incorporating feedback into 
semi-analytical and numerical models of galaxy formation. 

Fortunately it is possible to obtain some quantitative infor- 
mation about the impact of galactic winds without a full un- 
derstanding of the abstruse physics responsible for their mor- 
phology and kinematics. Galaxies which host winds power- 
ful enough to overcome the gravitational binding energy of 
their halos will vent some of their metals into the intergalac- 
tic medium. Hence, one way of evaluating the importance of 
galactic winds is to look for their chemical imprint on galax- 
ies. However, low metallicity is not necessarily a hallmark of 
wind activity. The metallicity of a galaxy is expected to de- 
pend strongly on its evolutionary state, namely how much of 
its gas has been turned into stars. To detect metal depletion 
it is therefore necessary to make some assumptions about the 
expected level of chemical enrichment based on a galaxy's 
star and gas content. Of course other mechanisms besides 
winds could make a galaxy appear to be metal-depleted — for 
example, the inflow of pristine gas, or the return of compara- 
tively unenriched material from evolved low mass stars. But 
these scenarios can potentially be distinguished by examin- 
ing the dependence of metal depletion on dynamical mass. To 
quantify the impact of feedback on the local galaxy popula- 
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tion we therefore compare the observed metallicities of galax- 
ies spanning a wide range in total mass to the predictions of 
simple chemical evolution models. 

Interest in the relationship between mass and metallicity 
da tes back severa l decad es, beginning with the seminal work 
of Leaueux et al] d!979h . However, because of the difficulty 
of obtaining masses, luminosity was generally adopted as a 
surrogate. A correlation b etween m etallicity and blu e lumi- 
nosity was demonstrated by Garnett & Shields ( 1987) and ex- 
tende d by various authors (e.g Skillman, Kennicutt, & Hodge 
| 1989[lBrodie & Huchrall99lUZaritskv. Kennicutt. & Huchral 
1994) to include a range of Hubble types and to span over 
11 magnitudes in luminosity and 2 dex in metallicity. The 
interpretation of this striking trend has, however, remained a 
matter of some debate owing to the difficulty of transforming 
from observables to the ingredients needed for simple chemi- 
cal evolution models (stellar mass, gas mass, and metallicity). 

In recent years, the developme nt of more sop histicated 
models for stellar popula tions (e.g. Bruzual & Chariot j^QQ^L 
and gaseous nebulae (e.g. Ferland 1996; Chariot & Longhetti 
120011) has resulted in major advances in our ability to derive 
physical properties from observables. In this work we de- 
rive st ellar mass-to-light (M/L) ratios according to the meth- 
ods of Kauffma nn et alJ fe003al) which rely on spectroscopic 
line indices to help circumvent the classical age-metallicity- 
reddening degene racy issues. We mea sure metallicity using 
the formalism of Chariot et al. (2004) which makes use of 
all of the strong optical nebular lines in our bandpass (as 
opposed to more traditional methods which rely on a single 
line ratio.) We couple these improved techniques with the 
enormous statistical power provided by the Sloan Digital Sky 
Survey (SDSS). We present the mass-metallicity relation of 
~ 53,000 star-forming galaxies at z ~ 0.1 as a new bench- 
mark for successful models of galaxy evolution. 

We begin with a brief description of the SDSS data prod- 
ucts and our data processing techniques in |2] We outline our 
method for measuring metallicity and stellar mass in 50 In |4] 
we compare the luminosity-metallicity relation of the SDSS 
data with previous determinations. We present the mass- 
metallicity relationship in ^5] and consider its origin in |6] 
We explore sources of systematic error in J7]and conclude in 
|8] Where appropriate we adopt a cosmology of fijyj = 0.3, 
Oa = 0.7, and H = 70 km s" 1 Mpc" 1 . 

2. THE SDSS DATA 

The data analyzed in this study are drawn from the Sloan 
Digital Sky Survey (SDSS). The survey goals are to obtain 
photometry of a quarter of the sky and spectra of nearly one 
millio n objects. Imaging is obtained in the u, g, r, i, z bands 
jFukueita et aljfl996t ISmith et alJl2 002l) w ith a special pur- 
pose drift scan camera (Gun n et al J 119 98) mounted on the 
SDSS 2.5 meter telescope at Apac he Point O bservatory. The 
imaging data are photometrica lly (Hogg et al] l2001l) and as- 
trometrically (Pieret al. 2003) calibrated, and used to select 
stars, galaxies, and quasars for follow-up fiber spectroscopy. 
Spectroscopic fibers are assigned to objects on the sky using 
an efficient tiling algorithm designed to optimize complete- 
ness (|Blanton et al. 2003a). The details of the survey strategy 
can be found in York etall fcOOOl) and an overview of the data 
pipeli nes and products is pro vided in the Early Data Release 
paper ( Stoughto n et all20 02). 

Our parent sample for this study is composed of 211,265 
objects which have been spectroscopically confirmed as 
galaxies and have data publicly available in the SDSS Data 



Release 2 (DR2; Abazai ian et alJl2004l) . These galaxies are 
part of the SDSS 'main' galaxy sam ple used for large scale 
structure studies ( Straus set alJ 120021) and have Petrosian r 
magnitudes in the range 14.5 < r < 17.77 after correction for 
foreground galactic extinction using the reddening maps of 
Schlegel. Finkbeiner. & Davis ( 1998). Their redshift distri- 
bution extends from ~ 0.005 to 0.30, with a median z of 0. 10. 
From this sample we select a subset of 53,4 00 st ar-forming 
galaxies for nebular analysis, as discussed in 32.21 

We determine stellar mass-to-light ratios and nebular metal- 
licities for our sample galaxies from the SDSS spectra. 
The spectra are obtained with two 320-fiber spectrographs 
mounted on the SDSS 2.5-meter telescope. Fibers 3" in di- 
ameter are manually plugged into custom-drilled aluminum 
plates mounted at the focal plane of the telescope. The spectra 
are exposed for 45 minutes or until a fiducial signal-to-noise 
(S/N) is reached. The median S/N per pixel for galaxies in the 
main sample is ~ 14. The spectra are processed by an auto- 
mated pipeline (Schlegel et al., in prep.) which flux and wave- 
length calibrates the data from 3800 to 9200 A. The instru- 
mental resolution is R = A/<5 A = 1 850 - 2200 (FWHM- 2 .4 A 
at 5000 A). 

The wavelength coverage, resolution, S/N, and general high 
quality of the SDSS spectra make them extremely well suited 
for the derivation of nebular abundances. However the spec- 
trophotometric calibration and the small size of the fiber aper- 
ture relative to the target galaxies pose some concern. The 
SDSS spectrographs do not employ an atmospheric dispersion 
corrector and the spectra are frequently acquired under non- 
photometric conditions. The Survey has nevertheless been 
able to obtain a remarkable level of spectrophotometric pre- 
cision by the simple practice of observing multiple standard 
stars simultaneously with the science targets. (The artifice in 
this case is that the 'standards' are not classical spectrophoto- 
metric standards, but are halo F-subdw arfs that are calibrated 
to stellar models - see Abazaii an et al.l J2004I) for details.) To 
quantify the quality of the spectrophotometry we have com- 
pared magnitudes synthesized from the spectra with SDSS 
photometry obtained with an aperture matched to the fiber 
size. The la error in the synthetic colors is 5% in g—r and 3% 
in r-i (X g ~ 4700 A; X r ~ 6200 A; A, - 7500 A). At the bluest 
wavelengths (~ 3800 A) we estimate the error to be ~ 12% 
based on repeat observations. There is also a systematic er- 
ror in the sense that the spectra are bluer than the imaging by 
^2% in the g-band, but it is unclear at present whether this 
represents an error in the absolute calibration of the photome- 
try or the spectroscopy. While the random errors in the spec- 
trophotometry are offset by the large size of the data set, the 
accuracy of our measured nebular abundances relative to true 
global abundances is compromised to some degree by the size 
of the fiber aperture. At the median redshift of the SDSS main 
galaxy sample (z ~ 0.1) the target galaxies subtend several 
arcseconds on the sky. The fraction of galaxy light falling in 
the 3" fiber aperture is typically about 1/3, with variations in 
redshift and morphological type causing the fraction to range 
from 0.05 - 0.6. The impact of this 'aperture bias' is discussed 
ind 

2.1. Emission Line Measurement 

The emphasis of this study is on the nebular emission line 
spectra of galaxies. However, we can not disregard the fact 
that the optical spectra of galaxies are very rich in stellar ab- 
sorption features, which can complicate the measurement of 
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nebular emission lines. In order to maintain speed and flexi- 
bility, the SDSS spectroscopic pipeline performs a very sim- 
ple estimate of the stellar continuum using a sliding median. 
While this is generally adequate for strong emission lines, a 
more sophisticated treatment of the continuum is required to 
recover weak features and to properly account for the stel- 
lar Balmer absorption which can reach equivalent widths of 
5 A in some galaxies. To address this need, we have de- 
signed a special-purpose code optimized for use with SDSS 
galaxy spectra which fits a stellar population model to the 
continuum. We adopt the basic assumption that any galaxy 
star formation history can be approximated as a sum of dis- 
crete bursts. Our library of template spectra is composed of 
single stellar population model s generated us i ng the new pop- 
ulation synthesis code of Bruzual & Chariot (2003, hereafter 
BC03). The BC03 models inc orporate an empirical spectral 
library (Le Bor gne et alJ 12003) with a wavelength coverage 
(3200 - 9300 A) and spectral resolution (~ 3 A) which is 
well matched to that of the SDSS data. Our templates in- 
clude models often different ages (0.005, 0.025, 0.1, 0.2, 0.6, 
0.9, 1.4, 2.5, 5, 10 Gyr) and three metallicities (1/5 Z Q , Z Q , 
and 2.5 Z Q ). For each galaxy we transform the templates to 
the appropriate redshift and velocity dispersion and resample 
them to match the data. To construct the best fitting model 
we perform a non-negative least squares fit with dust atten- 
uation modeled as an additional free parameter. In practice, 
our ability to simultaneously recover ages and metallicities is 
strongly limited by the signal-to-noise of the data. Hence we 
model galaxies as single metallicity populations and select the 
metallicity which yields the minimum \ 2 - (While this is not 
particularly physical, in practice it is not a bad assumption 
since the integrated light of a galaxy tends to be dominated 
by the light of its most recent stellar generation.) The details 
of the template fitting code will be presented in Tremonti et 
al. (in prep.). Figure^shows some typical spectra and their 
continuum models. 

After subtracting the best-fitting stellar population model 
of the continuum, we remove any remaining residuals (usu- 
ally of order a few percent) with a sliding 200 pixel me- 
dian, and fit the nebular emission lines. Since we are 
interested in recovering very weak nebular features, we 
adopt a special strategy: we fit all the emission lines 
with Gaussians simultaneously, requiring that all of the 
Balmer lines (H5, H7, H(3, and Ha) have the same line 
width and velocity offset, and likewise for the forbidden 
lines ([O II] AA3726,3729, [O III] AA4959,5007, [N II] 
AA6548,6584, [S II] AA6717,6731). We are careful to take 
into account the wavelength-dependent instrumental resolu- 
tion of each fiber, which is measured for each set of observa- 
tions by the SDSS spectroscopic pipeline from the arc lamp 
images. The virtue of constraining the line widths and veloc- 
ity offsets is that it minimizes the number of free parameters 
and effectively allows the stronger lines to be used to help 
constrain the weaker ones. Extensive by-eye inspection sug- 
gests that our continuum and line fitting methods work well. 

2.2. The Galaxy Sample 

From our original sample of ~ 2 1 1 , 000 galaxies, we select 
a subsample of star-forming galaxies for nebular analysis. We 
start by imposing a redshift cut of 0.005 < z < 0.25, and re- 
quiring that > 10% of the total galaxy light be observed by the 
fiber. Because stellar mass and metallicity are derived using 
Bayesian techniques (see 50, it would be possible to select a 
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FIG. 1. — Images and spectra of some typical galaxies in our sample at a 
range of redshifts. The red circle on the images denotes the fiber aperture. 
The continuum model for each spectrum is shown in red and offset slightly 
for cla rity. The continuum model is constructed from the Br uzual & Chariot 
1 2003 ) stellar population synthesis models as described in i|2.1l Both ob- 
served and model spectra have been smoothed by 5 pixels. 



sample of objects with accurate masses and metallicities us- 
ing the likelihood distributions. However, for clarity we have 
opted to define our sample based on observables (line fluxes, 
line indices) and later make some weak cuts based on the de- 
rived parameters. We require galaxies included in our star- 
forming sample to have lines of H/3, Ha, and [N II] A6584 
detected at greater than 5er. We do not explicitly constrain the 
other nebular lines that we make use of ([O II] AA3726, 3729, 
[O III] A5007, He I A5876, [O I] A6300, [S II] AA6717,6731) 
because our Bayesian analysis takes into account the errors 
in the line fluxes. We note that [O II] is not measured at 
Z < 0.03 due to the blue wavelength cut-off of the spectro- 
graphs. This affects 5% of our sample but has a negligible 
impact on our nebular analysis due to the use of multiple 
lines (see ^3. II . We impose the requirement that the param- 
eters needed for mass determination have small errors (see 
33. 2> . These criteria are: <r(m z ) < 0.15 mag, (t(HSa) < 2.5 A, 
and <t[D„(4000)] < 0.1. The combination of these constraints 
trims the sample to ~ 82, 000 galaxies, with the limiting factor 
typically being the strength of the H(3 line. 

We remove galaxies harboring an Active Galactic Nu- 
cleus (AGN) from our sample using the traditional 
line diagnostic diagram fN II]/Hq versus [O I II]/H/3 
( Baldwin. Phillips. & Terlevichll 19811 IVeilleux "& Osterbrock 
1987). The division between star- forming galaxies and AGN 
has been calibrated theoretically by Kewlev e t alj(|20 01[). We 
adopt th e slightly modified formula used bvlKauffman n et alJ 
(2003c) which provides a more conservative approach to se- 
lecting star-forming galaxies (see their Fig. 1). We apply 
this criterion where we detect [O III] with 3 a sigma signif- 
icance (and Hp, Ha, and [N II] at 5(7, as specified earlier.) 
To avoid biasing ourselves against high metallicity galaxies 
which have intrinsically weak [O III], we also include objects 
with [O III] < 3cr and log([N II]/Ha) < -0.4 in our star form- 
ing sample. These galaxies comprise ~ 16% of our final sam- 
ple. The fraction of emission line galaxies that we reject be- 
cause of AGN contamination is ^33%. 
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Finally, we require the la errors in stellar mass and metal- 
licity derived from the likelihood distributions to be less than 
0.2 dex. This cut eliminates 3% of the remaining galaxies, 
leaving us with a sample of 53,400 galaxies. Some example 
galaxies and their spectra are shown in Figure ^ The red- 
shift, luminosity, color, and morphology of our star-forming 
sample is compared with that of the DR2 main galaxy sam- 
ple in Figure [2] We use the concentration index, C, as a 
rough proxy for galaxy morphology. It is defined as the ra- 
tio of the radii enclosing 90% and 50% of the Petrosia n r- 
band galaxy light (C = R90/R50). Shimasaku et al. (200T|) and 
Strate vaet al J d200 1 ) have shown that there is a good corre- 
spondence between C and Hubble type, with C ~ 2.6 marking 
the boundary between early- and late-types. As expected, our 
emission line sample has bluer colors and more late-type mor- 
phologies than the parent sample. The redshift distribution is 
similar to that of the main survey, and peaks around z = 0.08. 




FIG. 2. — Properties of the SDSS main sample and our star-forming galaxy 
sample (gray shaded histograms). The upper left panel shows the redshift 
distribution. The upper right panel shows absolute r-band Petrosian magni- 
tudes. The lower left panel shows the distribution of gal axy (g — i) color . The 
colors have been ^-corrected to z = 0. 1 following Blanton et al. 1 2003b ). The 
lower right panel shows the concentration index, C = Rgo/Rso- Concentra- 
tion correlates loosely with Hubble type, with C ~ 2.6 marking the boundary 
between early- and late-types. 



3. THE PHYSICAL PROPERTIES OF GALAXIES 

3.1. Measuring Metallicity 

We derive oxygen abundances for our SDSS galaxies from 
the optical nebular emission lines. There are a number of ad- 
vantages of measuring metallicity 12 from the nebular lines 
rather than the stellar absorption features via Lick indices 
rWorthevlll99l . First, the S/N in the emission lines can 
greatly exceed that in the continuum, with the added advan- 
tage that many of the optically faintest galaxies (dwarfs) ex- 
hibit the highest emission line equivalent widths. Secondly, 
nebular abundance determination is free of the uncertainties 
due to age and a-element enhancement which plague the in- 
terpretation of absorption line indices. Thirdly, nebular metal- 
licities are easier to interpret in the context of chemical evolu- 
tion models because they reflect the present-day metal abun- 
dance rather than the luminosity-weighted average of previ- 



ous stellar generations. The disadvantage is that our analysis 
is limited to galaxies with on-going star formation. 

The strong optical nebular lines of elements other than H 
and He are produced by collisionally excited transitions. The 
strength of these lines is the product not only of the heavy el- 
ement abundance, but the temperature, density, and ionization 
state of the nebula, the attenuation by dust, and the depletion 
of metals on to dust grains. In practice, direct determination 
of element abundances from optical emission lines tends to be 
limited by the faintness of the [O III] A4363 line which is used 
to determine the electron temperature. When the electron 
temperature cannot be measured directly, oxygen abundances 
can be estimated using empirically or theoretically calibrated 
relations between metallicity and the relative fluxes of various 
strong optical emission li nes. 'Strong-line' abundanc e cali- 
bration was pioneered bv lAlloin et al.l (H979) and Pag el et all 
J1979I) . The latter introduced the R23 metallicity indicator 
(R23 = ([O II] A3727 + [O III] AA4959,5007) / H/3) which has 
enjoyed widespread use. Various other line ratios have been 
calibrated as metallicity indicators, including [N II] A6584/ 
Ha, [O III] A5007/[N II] A6584, [N II] A6584/[0 II] A3727, 
and ([S II] AA6717,6731 + [S HTj AA90069,9532)/H/3. Sys- 
tematic differences among t he various strong-line c alibrations 
can exceed 0.2 dex. (See iKewlev & D onita 2002 for a re- 
view). 

Here we measure metallicities in a slightl y more refined 
way, using the approach outlined by Cha riot et alJ (2004). 
This consists of estimating metallicity statistically, based on 
simultaneous fits of all the most prominent emission lines 
([O II], H/3, [O III], He I, [O I], Ha, [N II], [S II]) with a model 
design ed for the inte rpretation of inte grated galaxy s pectra 
(Chariot & Lon ghettliool . The lCharlot & Longhettil (120011) 
model is based on a combination of the Bruzual & Chariot 
i 19931) and Ferland (1996, version C90.04) population syn- 
thesis and photoionization codes. In this model, the contri- 
butions to the nebular emission by H II regions and diffuse 
ionized gas are combined and described in terms of an effec- 
tive (i.e. galaxy-averaged) metallicity, ionization parameter, 
dust attenuation at 5500 A, and dust-to-metal ratio. The de- 
pletion of heavy elements onto dust grains and the absorption 
of ionizing pho tons by dust are in cluded in a self-consistent 
way. Chariot & Longhetti (2001) calibrated the emission- 
line properties of their model using the observed [O III]/H/3, 
[O II]/[0 III], [S II]/Ha and [N II]/[S II] ratios of a representa- 
tive sample of 92 nearby spiral, irregular, starburst, and H II 
galaxies. The model also a ccounts well for the prop erties of 
these lines in our sample (Brinchmann et al. 2004, Fig. 5). 
We calculate the likelihood distribution of the metallicity of 
each galaxy in our sample, based on comparisons with a large 
library of models (~ 2 x 10 5 ) corresponding to different as- 
sumptions about the effective gas parameters. We adopt the 
median of this distribution as our best estimate of the galaxy 
metallicity; the width of the likelihood distribution provides 
a measure of the error. The median ler error for our sam- 
ple of star-forming galaxies is 0.03 dex. We note that this 

12 The word 'metallicity' has been used variously in the literature. In stel- 
lar studies metallicity usually refers to the iron abundance. Oxygen has been 
adopted as the canonical 'metal' for ISM studies because it is the most abun- 
dant, it is only weakly depleted on to dust grains, and it displays strong lines 
in the optical. Hereafter we use metallicity to denote the gas-phase oxygen 
abundance measured in units of 12 + log(0/H), where the ratio O/H is the 
abundance by number of oxyg en relative to hydrogen. Solar metallicity in 
these units is 8.69 1 Allende Prieto. Lambert. & Apslund 2001 ) and LMC and 
SMC metallicity are 8.4 and 8.0. respectively lGarnetlfl9 99). 
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does not include systematic error. Because our metallicities 
are discreetly sampled, we have added small random offsets 
for display purposes. These have a Gaussian distribution with 
a = 0.02 dex. 

Our models implicitly assume that the nebular gas 
in a galaxy can be characterized by a set of galaxy- 
averaged physical properties. An important issue is 
whether truly representative abundances can be derived 
for spatially integrated galaxy spectra in the presence 
of radial variations in temperature, ionization, metallic- 
ity, and dust extinction ty ila-Costas & Edmunds I1992J; 
Zarits kv. Kennicutt & Huchral 1 19941 Uvlartin & Rovl 1 19941 
vanZeeetal. 1998). Kobulnickv. Kennicutt. & Pizagno 
( 1999) addressed this issue by comparing H II region abun- 
dances at particular disk radii with abundances measured from 
integrated spectra. They found excellent agreement (±0.05 
dex) when metallicities were based on R23. However, the 
study is limited in one respect: they use weighted averages 
of the H II region spectra as proxies for true integrated spec- 
tra, thereby ignoring the contribution of the diffuse ionized 
gas which ca n be re sponsible for 25 - 70% of t he Balmer 
line emissi on ( Martin|1997|:|Zurita 1 R qzas 1 & Beckmanl2000t 
Thilker. Walterbos. Braun. & Hoopes 2002). Hence further 
study is needed to quantify the systematics associated with 
true integrated abundances. 

To place our metallicity measurements in the context of pre- 
vious work, we show the relation between our derived metal- 
licities and the line ratio R23 = ([O II] + O III)/H/3 in Figure^ 
For comparison we show t he theoretical rela tion between R23 
and metallicity derived by McGaush ( 1991 ) 13 , the empirical 
relation defined by Edmunds & Page] ( 1984 ), and the semi - 
empirical relation of Zaritskv, Kennicutt, & Huchra ( 1994), 
itself an average of 3 previous calibrations ( Donita & Evans 
1986; Edmunds & Pagel 11984 iMcCall. Rvbski. & Shields! 
1985). These analytic R23-metallicity relations roughly 
bracket the range of metallicities that we derive, showing 
that our measurements are in line with previous strong-line 
calibrations. However, recent work suggests that strong-line 
methods may overestimate ox ygen abundances systematically 
by as much as a factor of two ( Kennicutt. Bresolin. & Garnett 
2003). 

The strength of our Bayesian metallicity estimates is that 
they make use of all of the available nebular lines rather than 
relying on a small subset of them, as is the case for R23. How- 
ever, full spectral modeling is not always possible or practical. 
Hence we provide an analytical fit to the R23-metallicity rela- 
tion shown in Figure|3] 

12 + log(0/H) = 9.185 -0.313x-0.264x 2 -0.32Lr\ (1) 

where x = logR23. This formula is valid for the upper branch 
of the double-valued R23 -abundance relation. The lower 
branch is not very well sampled by our data. We measure 
metallicities below 12+log(0/H)=8.5 for 940 galaxies. How- 
ever less than 1/3 of these can be plotted in Figure |3]because 
most are at such low redshift that [O II] A3727 is not in our 
spectroscopic bandpass. While galaxies lacking [O II] have 
larger errors associated with their derived abundances, they 
show trends between abundance and luminosity that are gen- 
erally consistent with the rest of the sample. 

3.2. Measuring Stellar Mass 

13 We use the analytic formula for the relation between R23 and metallicity 
reported in Kobulnickv. Kennicutt. & Pizaano 1 1999). 
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log(R 3 = ([OII] + [OIII])/H|3) 

FIG. 3. — Comparison of the relation between metallicity and the line 
ratio R23 = ([O II] + [O III] / H/3). In this and subsequent figures we rep- 
resent the SDSS data using a combination of a two dimensional histogram 
and plotted points. We histogram the data where more than 5 data points 
fall in an individual pixel, and plot it as individual points otherwise. The his- 
tograms have been square-root scaled for better visibility. Metallicities for the 
SDSS data have been derived using statistical met hods describ ed in the text. 
The blue line shows the theoretical calibration of McGaugh 1 1 99 1 ) (as re- 
ported in Kobulnickv. Kennicutt. & Pizaeno 1 1999)) for three representative 
values of [O III]/[Q II]. The green line shows the empirical calibration of 
Edmunds & Pasel 1 1984), and the red line shows the semi-empirical calibra- 
tion of Zaritskv. Kennicutt. & Huchra 1 1994), itself the average of 3 previous 
calibrations. 



The stellar mass of a galaxy cannot be inferred directly from 
its optical luminosity because the stellar mass-to-light ratio 
depends strongly on the galaxy's star formation history and 
metallicity. Optic al colors have been widely used to estimate 
M/L ratios (e.g. iBell & de Jongll2001t iBrinchmann & Ellisl 
2000). However, errors in the derived M/L are known to re- 
sult if galaxies have formed a substantial fraction (> 10%) of 
their stars in a recent burst. This caveat is particularly prob- 
lematic given our sample of activ ely star-forming g alaxies. 
We therefore adopt the method of iKauffmann et alT f2003a) 
for deriving stellar masses. This method relies on spectral in- 
dicators of the stellar age and the fraction of stars formed in 
recent bursts. As our present sample is actually an extension 
of Kauffmann's, we provide only a cursory description of our 
methods here. 

We use the z-band (A ~ 8900 A) magnitude to characterize 
the galaxy luminosity because it is less sensitive than the bluer 
bands to extinction and the age of the stellar population. How- 
ever, even in this near-infrared passband the M/L ratio is not a 
constant, but can vary by up to a factor of 1 when a full range 
of star formation histories is considered (IKauffma nn et alJ 
20033 see Fig. 4). Constraints on the star formation history 
are provided by the spectral indices D„(4000) and HSa, which 
measure the 4000 A break and the stellar Balmer absorption. 
The location of a galaxy in the D„(4000)- HSa plane is insen- 
sitive to reddening, and only weakly dependent on metallicity; 
however, it is a powerful diagnostic of whether the galaxy has 
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been forming stars continuously or in bursts over the past 1 
- 2 Gyr. We assign stellar M/L ratios to our galaxies by us- 
ing a Bayesian analysis to associate the observed D„(4000) 
and HSa values with a model drawn from a large library of 
Monte Carlo realizations of galaxies with different star for- 
mation histories and metallicities. A comparison with broad- 
band photometry then yields estimates of the dust attenuation 
and the stellar mass. Errors in our masses are estimated from 
the width of the likelihood distribution. The median la er- 
ror for our samp le of star -forming galaxies is 0.09 dex. Our 
masses assume a Krourja (2001) Initial Mass Function (IMF). 

4. THE LUMINOSITY-METALLICITY RELATIONSHIP 

Because of the relative difficulty of measuring the 
stellar mass of galaxies, most previous work has fo- 
cused on the relationship between galaxy luminosity 
and metallici ty. Fo llowi ng the p ioneerin g work of 
iLeaueux etafl ( I19791) . ISkillman. Kennicutt. & Hodgd 
established a luminosity-metallicity rela- 
tion for irregul a r gal axies which was confirmed by 
iRicher & McCalll ill 9951) and e xte nded to spiral galaxies 
by iGarnett & Shields! J 19871). IVila-Costas & Edmunds! 
i 19921). IZaritskv. Kennicutt. & Huchral i 1994ft . and 
Kobulnickv & Zaritskv ( 1999). Recently attempts have 
been made to provide better comparison samples for high 
redshift by studyin g the luminosity-metallicity re lation 
in st arburst nuclei llCoziol. Contini. Davoust. & Considerd 
119971) . and UV-select ed (ICon tini et alJ 120021) an d emission 
line-selected samples (Melbourne & Salzer 2002). 

It is instructive to compare the luminosity-metallicity re- 
lation of our sample with previous work, both to validate 
the consistency of our measurements and to elucidate the ef- 
fects of our sample selection. Although redder bandpasses 
are less sensitive to the effects of dust obscuration and recent 
starbursts, historical precedent mandates that the luminosity- 
metallicity relationship be presented in terms of absolute B 
magnitude. Fortunately, the SDSS g band is fairly similar to 
the Johnson B ban d. We adopt the transformation given in 
ISmith etal] (1200 2) and fc-co rrect the magnitudes to z = us- 
ing the kcorrect code of Blanton et al. (2003b). We cor- 
rect our m easured blue luminosities for galact ic foreground 
extinction ( Schlegel, Finkbeiner, & Davis 1998) and for incli- 
nation dependent intrinsic attenuation ( Tully et al. 1998J). The 
latter implies a correction to a face-on orientation, but does 
not account for the intrinsic attenuation in a face-on system. 
The median correction is A%° = 0.3 mag. 

Figure [4] shows the luminosity-metallicity relation of the 
SDSS galaxies and various samples from the literature. We 
have corrected the published B-band luminosities to our 
adopted value of Hq where appropriate, but we have made 
no attempt to homogenize the metallicities to account for the 
different calibrations used. In comparing SDSS data with 
other samples, it is important to consider the possible aper- 
ture bias. We expect our metallicity measurements to slightly 
overestimate the true global abundances, as discussed further 
in 50 However, in view of the different metallicity calibra- 
tions used, the magnitude of this offset is probably not the 
dominant systematic effect. The SDSS data show generally 
good agreement with other local samples of galaxies, albeit 
with a paucity of dwarfs. 

In keeping with previous work, we fit the luminosity- 
metallicity relation using a linear least squares technique. The 
measurement errors in luminosity (±0.01 dex) and metallicity 
(±0.1 dex) are small compared to the observed scatter. The 



traditional method of estimating the underlying functional re- 
lation in this case is to use the linear bisector — the line 
which bisect s the ordinary least-squares regression of Xvs. Y 
and Yvs. X (llsobe. Feigelson. Akri tas. & Babu 1990). The 
luminosity-metallicity relation for our galaxies fitted in this 
manner is 

12 + log(O/H) = -0.185(±0.001)M B + 5.238(±0.018). (2) 

Our measured slope is interm ediate between that found by 
Kobulnickv & Zaritskv ( 1999 ) for local irregula r s and spi- 
rals and that measured bv iMelbourne & Salzerl ( 120021) for 
a subset of galaxies from the KPNO International Spec- 
troscopic Survey (KISS). In inter-comparing data sets it is 
important to note that different authors (including the two 
mentioned above) have adopted different regression meth- 
ods. Different regression methods calculate intrinsically dif- 
ferent p roperties of the data and are not directly compara- 
ble (see llsobe. Feigelson. Akritas. & Babul (119901) for a full 
discussion). Figure |4] therefore elucidates the differences in 
the measured luminosity-metallicity relation which are the 
result of sample selection, metallicity calibration, and data 
characterization (fitting). All of these factors need to be care- 
fully considered when looking for evidence of evolution in the 
luminosity-metallicity relation. 
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FIG. 4. — The luminosity-metallicity relation for SDSS galaxies and var- 
ious galaxy samples drawn from the literature (see legend). All of the B 
magnitudes have been corrected to Hq = 70 km s Mpc~' , but we have made 
no attempt to homogenize the metallicity measurements. The red line rep- 
resents the linear least squares bisector fit to the SDSS data. The inset plot 
shows the residuals of the fit. 

At fixed metallicity, the SDSS galaxies span more than 4 
magnitudes in blue luminosity. The brightest galaxies are 
consistent wi t h the l uminosity-metallicity relation found by 
IContini et alJ (|2002) for a sample of UV-selected galaxies 
(pink line in Fig.^i. Contini et al. postulate that these galaxies 
are offset from the luminosity-metallicity relation of normal 
spirals due to the low M/L ratio of their newly formed stellar 
populations. This suggests a physical origin for the scatter, an 
idea that we explore further in Figure|5] To reduce sources of 
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measurement error, we switch to using the luminosity in the 
native SDSS passbands, ^-corrected to the median redshift, 
Z ~ 0.1. To quantify the influence of dust, we correct the 
galaxy luminosi ties for intrinsi c attenuation using the atten- 
uation curve of ICharlot & Falll ll200 0) and assuming that the 
stars experience 1/3 of the attenuation measured in the neb- 
ular gas. The nebular attenuation is determined simultane- 
ously with the metallicity assumin g a metallic ity-dependent 
Case B HafHf3 ratio JBrinchmann et al.ll20"oH ICharlot et alJ 
2004). We assume that the correction for intrinsic attenuation 
automatically accounts for any inclination-dependent effects. 

In Figure [5] we systematically examine the impact of dust 
and M/L variations on the luminosity-metallicity relation by 
adopting different measures of galaxy luminosity: 1) the ab- 
solute g-band magnit ude corrected for in clination-dependent 
attenuation following Tu llv et al.l (119981) : 2) the absolute g- 
band magnitude corrected for intrinsic attenuation, as de- 
scribed above; 3) and the absolute z-band magnitude cor- 
rected for intrinsic attenuation. In panels 1-3 of Figure |5] 
we indicate the distribution of metallicity at a given lumi- 
nosity by displaying the contours which enclose 68 and 95% 
of the data in bins of 0.4 mag. The contours provide a non- 
parametric description of the distribution which is unbiased 
as long as the errors in luminosity are small relative to our 
adopted bin-size. For comparison, we also show the tradi- 
tional least-squares linear bisector fit to the data in panel 1 
(12+log(0/H) = -0.186(±0.001)M ? + 5.195(±0.018)). Be- 
cause we do not know a priori the true functional form of 
the luminosity-metallicity relation, we focus on the contours. 
Comparison of the first two panels of Figure|5]shows that cor- 
recting the luminosity for attenuation reduces the scatter and 
flattens the luminosity-metallicity relation at high mass. This 
trend is even more pronounced when the extinction corrected 
z-band magnitude is used (panel 3). Because the z-band is 
less sensitive to dust and recent starbursts, the range of M/L 
ratios is smaller, and the scatter is reduced by ~ 20% com- 
pared to the uncorrected g-band. However, even in the z-band, 
M/L ratios can vary by factors of a few. This effect is illus- 
trated in panel 4 where we plot the median z-band luminosity- 
metallicity r elation for galaxie s in four bins of D„(4000). As 
discussed in Kauffmann et al. (2003a), D„(4000) is a good 
measure of the mean stellar age of the population. Our inter- 
pretation of panel 4 is that at fixed metallicity, galaxies with 
lower D„(4000) are shifted to brighter magnitudes because of 
the lower M/L ratios of their young stellar populations. This 
confirms our intuition that the underlying physical correlation 
is between stellar mass and metallicity. 

5. THE MASS-METALLICITY RELATIONSHIP 

With our new prescriptions for measuring stellar mass and 
gas-phase metallicity it is now possible to examine the mass- 
metallicity relationship of our sample of SDSS star-forming 
galaxies. Figure |6] shows that a striking correlation is ob- 
served, extending over 3 decades in stellar mass and a fac- 
tor of 10 in metallicity. The correlation is roughly linear from 
10 8 ' 5 M Q to 10 M Q after which a gradual flattening occurs. 
Most remarkable of all is the tightness of the correlation: the 
lcr spread of the data about the median is ±0.10 dex, with 
only a handful of extreme outliers present. The relationship is 
well fitted by a polynomial of the form: 

12+log(0/H) = -1 .492 + 1 .847(logM«) - 0.08026(logM*) 2 

(3) 

where M* represents the stellar mass in units of solar masses. 
This equation is valid over the range 8.5 < log M» < 11.5. 
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FIG. 5. — The luminosity-metallicity relation of SDSS galaxies in the g and 
z-bands. In the first panel we have corrected M g to face-on orientation, but 
we have not corrected for internal attenuation. In panels 2-4 we correct M g 
and M- for internal attenuation, assuming that the stars experience 1/3 of the 
reddening measured in the gas. The solid black contours in panels 1-3 enclose 
68 and 95% of data with statistics computed in bins of 0.4 mag in luminosity. 
The median half-width of the distribution is listed in the lower right corner. 
For comparison, the dashed line in the first panel shows the least-squares 
linear bisector fit to the data. The fourth panel shows the median z-band 
luminosity-metallicity relation for galaxies in four bins of D„(4000): from 
bottom to top, 1.0 - 1.2, 1.2 - 1.3, 1.3 - 1.4, 1.4 - 1.8. Data for the contours in 
panels 1 and 3 are given in Tables 1 and 2. 
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FIG. 6. — The relation between stellar mass, in units of solar masses, and 
gas-phase oxygen abundance for ~53,400 star-forming galaxies in the SDSS. 
The large black points represent the median in bins of 0. 1 dex in mass which 
include at least 100 data points. The solid lines are the contours which enclose 
68% and 95% of the data. The red line shows a polynomial fit to the data. 
The inset plot shows the residuals of the fit. Data for the contours are given 
in Table 3. 
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and luminosity-metallicity relations is the more pronounced 
turnover seen at high metallicitywhen mass is used as the in- 
dependent variable. In Figure |5] we demonstrated that both 
dust and M/L variations act to smear out the turnover when 
luminosity is used in place of mass. Accordingly, the scat- 
ter in the uncorrected g-band luminosity-metallicity relation 
is ~50% higher than the scatter in the mass-metallicity rela- 
tion. 

We emphasize that the turnover in the mass-metallicity re- 
lation is not an artifact of our metallicity calibration. As 
shown in Figure[3]our metallicity estimates are consistent with 
other strong-line calibrations. If we adopt the R73 me tallic - 
ity calibration of Zaritsk vTKennicutt. & Huchral (ll994l) . the 
shape of the mass-metallicity relation is well r eproduced, 
but th e scatter increases by ~ 15%. Use of the McGaugh 
( 119911) R23 metallicity calibration results in a more pro- 
nounced turnover, due to the flatter relation between R23 and 
log(0/H). 

The shape of the mass-metallicity relation is also robust 
with respect to our choice of S/N cuts. The S/N in the H/3 
line is generally the parameter that drives our galaxy selec- 
tion. The main concern is that our requirement of (S/N)h/3 > 5 
may bias us toward massive galaxies with abnormally high 
gas fractions. However, at fixed equivalent width (EW) the 
distribution of S/N values is broad. Galaxies with H(3 EW 
< 5 A have a median (S/N)//^ of 1 1, and more than 10% of 
them have (S/N)/^ > 20. Thus our cuts on S/N are not im- 
posing stringent cuts on EW or gas fraction. We have also 
done the straightforward experiment of altering the S/N cuts 
on all our lines and re-fitting the mass-metallicity relation: its 
detailed shape is insensitive to factor of 2 changes in our S/N 
cuts. 

The most remarkable feature of the stellar mass-metallicity 
relation is the tightness of the correlation - the la spread of 
the data about the median ranges from 0.20 dex at low mass to 
0.07 dex at high mass. For comparison, the median errors in 
mass and metallicity are 0.09 dex and 0.03 dex respectively. 
Adding these appropriately, we find that roughly half of the 
spread in the mass-metallicity relation can be attributed to 
observational error. With our large statistical sample we can 
test whether there is a physical origin for the scatter, even if 
we have underestimated our measurement errors. In Figure 
we show the residuals from our fit to the mass-metallicity re- 
lation as a function of various galaxy physical properties. 

There is a clear tendency for the metallicity residuals to 
correlate with local surface mass density measured within 
the fiber aperture. Analogous trends were identified by 
Bell & de Jongl ( 120001) using resolved photometry of a sample 
of nearby disk galaxies. They demonstrated that mean stellar 
metallicity and stellar age correlate with local K-band surface 
density and concluded that surface density plays a strong role 
in shaping both the local and global star formation history of a 
galaxy. Kauff mann et alJ (2003b) showed that the star forma- 
tion histories of low mass galaxies are more fundamentally 
related to surface mass density than to stellar mass. These 
findings suggest that at fixed mass, galaxies with higher sur- 
face densities have transformed more of the available gas into 
stars, thereby raising the gas-phase metallicity. This interpre- 
tation is borne out by the fact that the residuals also correlate 
somewhat with galaxy color. But curiously, the same trend is 
not seen with Ha EW or galaxy morphology as measured by 
the concentration index. This leaves us with a somewhat con- 
tradictory picture since these quantities are expected to relate 



to a galaxy's gas content and general evolutionary state. 

The other pronounced correlation seen in Figure is be- 
tween metallicity and inclination, as measured by the ratio of 
the galaxy minor and major axes. At fixed mass, fully edge-on 
galaxies (b/a~ 0.2) are more metal poor than face-on galaxies 
by 0.19 dex. This trend is easily explained because more of 
the disk is seen in projection in our fiber aperture at higher 
inclinations. The integrated metallicity decreases as larger 
galaxy radii are probed because the outer regions of disks tend 
to be more metal-poor, as well as more gas rich and less ex- 
tincted. 
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FIG. 7. — Correlation of the residuals of the mass-metallicity relation with 
various galaxy physical parameters: A/*, the stellar mass in units of Mq; 
the local surface mass density measured within the fiber aperture in units of 
Mq kpc~ 2 ; C, the concentration parameter, a proxy for galaxy morphology; 
b/a, the ratio of the galaxy minor and major axes, a proxy for inclination; 
(g — i) color ^-corrected to z=0. 1 ; and Ha equivalent width in A. The colored 
lines show the median relation in 5 different mass bins. 



6. THE ORIGIN OF THE MASS-METALLICITY RELATION 

The existence of a tight (±0.1 dex) correlation between 
stellar mass and metallicity reflects the fundamental role that 
galaxy mass plays in galactic chemical evolution. However, 
it is not a priori clear whether this sequence is one of enrich- 
ment or of depletion. Simply put, if more massive galaxies 
form fractionally more stars in a Hubble time than their low- 
mass counterparts, then the observed mass-metallicity rela- 
tion represents a sequence in astration. However, if galaxies 
form similar fractions of stars, then the relation could imply 
that metals are selectively lost from galaxies with small po- 
tential wells via galactic winds. 

Both of these ideas have a long history in the astronomical 
literature. It is well known that there are systematic trends 
in the star fo rmation history of galaxie s along the Hubble se- 
quence (e.g. Roberts & Havneslll99l . Recent observations 
have confirmed t hat gas mass fractions decrease with increas- 
ing stellar mass dMcGaugh & de Blokll 1 997t iBell & de Jong l 
2000; Boselli, Gavazzi, Donas, & Scodeggio 2001), a trend 
that seems to suggest that low mass galaxies are un-enriched 
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rather than depleted. However, observations of starbursts 
have revealed the nearly ubiquitous presence of galactic winds 
jHeckmarJ 2002). while studies of X-ray bright clusters have 
demonstrated th e presence of copious metals in the intra - 
cluster medium ( Gi bson. Loewenstein. & Musho tzkvll 19971) . 
and absorption line studi es have revealed metals in the inter - 
galactic medium (Ellison, Sonaaila, Schave, & Pettini 2000). 
Given the existence of these contradictory pieces of informa- 
tion, anecdotal arguments are of little use: a direct test of the 
origin of the mass-metallicity relation is required. 

The effects of astration and mass loss on the relationship 
between mass and metallicity can be disentangled in princi- 
ple if information on the relative mass of the gas and stars 
is available. Simple closed box chemical evolution models 
predict that metallicity (Z) is a simple function of the stellar 
yield, y, and of the gas mass fraction, /j, gas : 

Z = yln(p- l as ) (4) 

Following IGarnettl (12002). we assume that the stellar yield 
is constant. Metallicity is then straightforwardly related to 
the gas mass fraction if the tenets of the simple model ap- 
ply - namely that there are no gas inflows or outflows. We 
can invert Equation 1 to define the 'effective yield', which 
can be computed from the observed metallicity and gas mass 
fraction. When the simple model applies, the effective yield 
will equal the true yield, independent of galaxy mas s. This 
straigh tforward observational test was performed bv IGarnettl 
( 2002) using a sample of 44 nearby spiral and irregular galax- 
ies. Garnett found the effective yield to be constant for galax- 
ies with rotational velocities in excess of ~ 150 km s _I and to 
decline by a factor of ~ 15 below this threshold. 

While our SDSS dataset is much larger than the nearby 
sample of Garnett, we do not have direct information about 
the H I and H2 gas content of our galaxies. However, we do 
have indirect information, by virtue of the fact that all of our 
sample galaxies are actively forming stars. To estimate the 
gas mass, we invoke another well k nown empirical corre lation 
- the Schmidt star formation law ( Schmida ll959t iKennicutil 
1998) which relates the star formation surface density to the 
gas surface density. 

For each of our galaxies we calculate the star formation rate 
(SFR) in the fiber ap erture from the atten uation-corrected Ha 
luminosity following Brinch mann et al.l ll2004 j ). We multiply 
our SFRs by a factor of 1.5 to conv ert from a|K roupa (2001) 
IMF to the Salpeter IMF used by iKennicuttl ( 119981) . Our 
SDSS galaxies have star formation surface densities which 
are within a factor of 10 of ^sfr = -3 M Q yr" 1 kpc -2 , exactly 
the range found by Kennicutt ( 1998) for the central regions of 
normal disk galaxies. We convert star formation surface den- 
sity to surface gas mass density, by inverting the com- 
posite Schmidt law of iKennicuH 11998). 



E SF r= 1.6 x 10" 



1M pc" 



1.4 



M Q yr 'kpc 2 . 



(5) 



(Note that the numerical coefficient has been adjusted to in- 
clude helium in Ti gas .) Combining our spectroscopically de- 
rived M/L ratio with a measurement of the z-band surface 
brightness in the fiber aperture, we compute S J/ar , the stel- 
lar surface mass density. The gas mass fraction is then /i gas = 

In Figure [8] we plot the effective yield of our SDSS star 
forming galaxies as a function of total baryonic (stellar + 
gas) mass. Baryonic mass is believed to correlate with dark 



mass, as evidenc ed by the existence of a baryonic ' Tully- 
Fisher ' relation llMcGaugh. Schombert. Bothun. & de Blokl 
120001 iBell & de Jondl2001l) . We are interested in the dark 
mass because departures from the 'closed box' model might 
be expected to correlate with the depth of the galaxy potential 
well. Because very few of our SDSS galaxies have masses 
below 10 85 Mr?), we augment our data s e t with measure - 
men ts from iLee. McCall. & Richer] ( 120031) . IGarnettl ( 120021) . 
and Pilvugin & Ferrini (2000), all of which use direct gas 
mass measurements. We note that the correspondence be- 
tween the SDSS data and the samples from the literature is 
very good in spite of the fact that the gas masses have been 
determined in very different ways. 
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FIG. 8. — Effective yield as a function of total baryonic mass (stellar + 
gas mass) for 53,400 star-forming galaxies in the SDSS. The large black 
points represent the median of the SDSS data in bins of 0.1 dex in mass 
which include at least 100 data points. The solid lines are the contours 
which enclose 68% an d 95% of the data. T he colored crosses are data from 
ILee. McCall. & Richer! 12001 IGarnenl 120021 and IPilvugin & FemnJ 1200(1 
Both the metallicities and the gas masses used to derive the effective yield 
have been computed differently in the SDSS data and the samples from the 
literature. The agreement nevertheless appears quite good. The pink line is 
the best fit to the combined dataset assuming the intrinsic functional from 
given by Equation 6. The dashed line indicates yo, the true yield if no metals 
are lost, derived from the fit to the data. The pink stars denote galaxies which 
have lost 50% and 90% of their metals. Data for the contours are given in 
Table 4. 

With the addition of the data from the literature it is clear 
that galaxies do not evolve as 'closed boxes'. Figure|8]shows 
that baryonic mass and effective yield are highly correlated, 
with the effective yield decreasing by a factor of ~ 10 from the 
most massive galaxies to dwarfs. The correlation is steep at 
low masses but begins to flatten around 10 9 5 M Q . In 5j8]we ar- 
gue that the relationship between effective yield and baryonic 
mass is the consequence of metal loss via galactic winds. We 
fit the combined dataset with a simple function predicated on 
this basis. We assume that the fraction of newly synthesized 
metals retained by a galaxy is proportional to the depth of 
the galaxy's potential well for low mass galaxies, but asymp- 
totes to unity for the most massive systems. The depth of the 
galaxy potential well scales approximately as V t 2 , where V c is 
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the galaxy circular velocity. We model the retained fraction as 
fret = V C 2 /(V C 2 +V$) where V is a constant. (f ret (V c < Vq) < 0.5, 
fret(Vc Vq) ~ 1 .) We ad opt M^vm oc V" r 3,5 fro m the baryonic 
Tully-Fisher relation of Bell & de Jong (2001). The effective 
yield and the baryonic mass may then be related follows: 



yeff = 



yo 



(l + (M /M baryon fV) 



(6) 



The free parameters of the fit are yo, the true yield (if no met- 
als are lost), and Mo, the mass at which a galaxy loses 1/2 of 
its metals. We perform the fit to a combination of the SDSS 
data (the median value of y e ff in bins of 0.1 in Mt, aryon ) and 
the data drawn from the literature, applying weights such that 
both datasets contribute equally to the fit. We find yo = 0.0104 
and Mo = 3.67 x 10 9 M Q . This characteristic mass corre- 
spo nds to a rotation speed of V r = 85 km s " 1 and , follow- 
ing Heckman. Lehnert. Strickland. & Armus (2000), an es- 
cape velocity of V esc ~ 3V C = 260 km s" 1 . Thus, the escape 
velocity of a galaxy that has lost half of its metals is just below 
the observed termin al velocity of starburst winds, which are 
300 - 900 km s _I (Heckman. Lehnert. Strickland. & Armus 
2000). 

The scatter in the correlation between effective yield and 
baryonic mass is ±0.15 dex, 50% larger than that measured 
for the mass-metallicity relation. However, much of the scat- 
ter in the SDSS data can probably be ascribed to our relatively 
uncertain measurement of the gas mass. We note that this does 
not effect the total mass greatly, as the median gas mass frac- 
tion of our sample is ~20%. It is interesting to note that galax- 
ies with direct gas mass measurements show comparable scat- 
ter, however the stellar mass determinations for these galaxies 
are correspondingly more uncertain. It is therefore premature 
to consider a physical origin for the scatter, although it is quite 
likely that the re may be one, for instance differences in galaxy 
environment (Skillman, Kennicutt, Shields, &Zaritskv 1996; 
iLee. McCall. & Richerl2003l) . 

7. SOURCES OF SYSTEMATIC ERROR 

Galaxies are know n to have radial grad ients in their phys- 
ical properties (e.g. IB ell & de Jon g| 120001) . Aperture effects 
are therefore a potentially serious concern because strong se- 
lection effects are inherent in magnitude-limited surveys such 
as the SDSS. For our sample, the median projected fiber size 
is ~4.6 kpc in diameter, with the range extending from 1 to 
12 kpc. Thus while we are not measuring 'nuclear' spectra, 
we are clearly not measuring integrated spectra either. The 
median fraction of galaxy light in the fiber aperture is 24%. 
Fortunately, because larger galaxies are brighter, they get se- 
lected out to larger distances and observed with larger pro- 
jected apertures. The net effect is that the fraction of galaxy 
light seen by the fibers is not a strong function of absolute 
magnitude, as shown in the top panel of Figure [9] Hence, 
while aperture effects are important in the absolute sense, they 
are unlikely to have a significant effect on trends with galaxy 
mass. 

Because of the magnitude-limited nature of the SDSS, it 
is not possible to examine the impact of aperture size sim- 
ply by comparing galaxies in narrow bins of absolute mag- 
nitude at different redshifts. Instead we divide our sam- 
ple into four broad bins of absolute magnitude and examine 
trends in galaxy properties as a function of relative redshift, 
z/Zmax, where Zmax is the redshift at which a given galaxy 
reaches the spectroscopic survey limits. The effect of aper- 
ture bias on galaxy M/L ratios was examined as a function 



of zjzmax bv lKauffmann et all d2003al see Fig. 18) and found 
to be ^0.1 dex. The strongest trends were observed for 
L* galaxies which presumably have both a well developed 
bulge and a disk. The effect of aperture bias on our mea- 
sured metallicity is shown in Figure |5] We find a change of 
at most -0. 1 1 dex in metallicity as galaxies move from one 
edge of the survey to the other, indicating that our metallic- 
ity determinations are only moderately affected by changes in 
the projected aperture size. This gradient is in good accord 
with our expectations based on the known radial m etallic- 
ity gradients in spiral galaxies jVila- Costas & Edmunds 1992; 
Zarits kv. Kennicutt & Huchrai 119941 iMartin & Rovl 1 19941 
Ivan Zee et alll998l) . 

While it is desirable to correct our metallicity measure- 
ments to reflect true global abundances, it is not possible to 
derive this correction directly from the SDSS data because 
even at the limits of the survey the fiber aperture only contains 
about 50% of a galaxy's light. Work is under way to obtain the 
observations necessary for such an analysis. In the meantime, 
we simply note that because the fibers generally cover the in- 
ner few kpc of galaxies, the abundances we derive are likely 
to be higher than true integrated abundances. Recent evidence 
also suggests that oxygen abundances derived from 'strong 
lines' methods like ours may overestimate the true abundance 
by as much as factor of two ( Kennicutt . Bresolin. & Garnettl 
2003). 



u 80 

D 

E 60 
g 

2 40 

a 

- 20 




-24 



-22 



-20 
M 



-18 



-16 



9.4 
9.2 

X 

O 9.0 



2 8.6 
8.4 
8.2 
9.4 

9.2 

s 90 

60 

3 8.8 
2 8.6 

8.4 
8.2 





-24 < M / < -22 




-22 <M 


<-21 












_ A log(Offl) = 


0.05 dex 

' ' 1 ' ' ' 1 ' ' ' 1 ' 


. A log(0/H) = 
' 1 ' ' ' 1 ' 


0.09 dex 


''I 




-21 <M <-20 




-20 <M 


<-18 










Hi ~ 


_ Alog(0/H) = 


0. 1 1 <lcx 


. Alog(0/H) = 


0.09 dex . 





0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

z / z z / z 

FIG. 9. — A test for aperture bias. In the top panel the fraction of r-band 
galaxy light observed by the fiber is plotted as a function of absolute mag- 
nitude. In the lower panel, metallicity is plotted in 4 broad bins of absolute 
magnitude as a function of relative redshift, z/zmax, where Zmax is the redshift 
at which a given galaxy reaches the survey limits. The large black points 
represent the median; the solid line is a linear fit to the points. The maxi- 
mum change in metallicity with z/zma (based on the fitted line) is given in 
the lower left of each panel. 



Another potential source of error is our indirect method of 
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deriving the gas mass from the H a luminosity via the global 
Schmidt law of Kennicutt ( 1998). While the impressive corre- 
lation between the surface mass density of gas and the surface 
density of star formation holds over 5 orders of magnitude, 
there is a fair amount of scatter evident in the na rrow regime 
applicable to the present study (Kennicutt 1998, see Fig. 6) . 
Systematic errors are a bigger concern. Won g & Blitz! 12002) 
found that application of a radially-dependent attenuation cor- 
rection to Ha produced a steeper Schmidt law, with a power 
law index of 1.7 (as compared to Kennicutt's 1.4). However, 
adopting the steeper Schmidt law in our analysis produces a 
relatively small change in the correlation between baryonic 
mass and effective yield: logyo decreases by 0.04 dex, and 
logMo decreases by 0.15 dex. 

Beyond the issue of the accuracy of our fj, gas measurement, 
there is the larger question of whether the central gas mass 
fraction is the relevant gas mass fraction. If T, gas /T, star varies 
strongly with radius, than our measurements are unlikely to 
be representative of the galaxy as a whole. These concerns 
are allayed to some degree by the generally strong correspon- 
dence s een between Ha and broad-band optical d isk scale 
lengths JHodge & Kennicutll983llKennicuttlll989l) . A more 
fundamental question, however, is how the products of star 
formation are distributed in galaxies - in other words, how 
well mixed is the interstellar medium? The mere existence 
of radial metallicity gradients in spirals suggests that mixing 
timescales are generally long. In this case, local abundances 
reflect local astration levels, and hence the relevant gas mass 
fraction is the local one. A final issue is the existence of cold 
gas external to the stellar disk in some galaxies. It is unclear 
whether this gas is relevant to our calculation, since it does 
not strongly participate in the star formation. However, if 
the outlying gas viscously mixes with gas in the inner disk it 
would affect the observed metallicity, and should be included 
in the gas budget. In any case, since we have no information 
on the extended gas in our galaxies, we do not consider it at 
present, except to note that its inclusion would raise our mea- 
sured value of the effective yield. For galaxies with baryonic 
masses of 10 7 M Q the gas mass would need to increase by a 
factor of ^ 1 1 to erase the signatures of metal loss. 

While the sources of systematic error cannot be ignored, we 
believe the magnitude of these errors to be small, on the order 
of a few tenths of a dex. In view of the moderately large dy- 
namic range of our results - a factor of ~ 10 decrease in both 
metallicity and effective yield - our general result, the prefer- 
ential loss of metals from galaxies with small potential wells, 
appears very robust. Future observations which spatially re- 
solve the distribution of gas and metals in galaxies will prove 
key to understanding many of the relevant systematic effects. 

8. SUMMARY & DISCUSSION 

We have coupled advanced techniques for deriving stellar 
masses and gas-phase metallicities from optical spectroscopy 
and photometry with the statistical power provided by the 
SDSS. For a sample of 53,400 star-forming galaxies at z ~ 
0.1, we find a tight correlation (±0.1 dex) between stellar 
mass and the gas-phase oxygen abundance, which extends 
over 3 orders of magnitude in stellar mass and a factor of 10 in 
oxygen abundance. We use indirect estimates of the gas mass 
fraction based on our measured Ha luminosity to estimate the 
effective yield. Simple closed box chemical evolution mod- 
els predict that the effective yield is constant. Combining our 
SDSS data with data from the literature, we find evidence that 
the effective yield decreases by a factor of ~ 10 from the most 



massive galaxies to dwarfs. 

The most straightforward interpretation of the correlation 
between baryonic mass and effective yield is the selective 
loss of metals from galaxies with shallow potentia l wells 
via ga lactic winds, an idea first introduced by iLarsonl 
( 1974). Evidence for galactic outflows has stead ily accumu- 
lated from observ ations at optical (|Lehnert & Hec kman 
119961 IHeckman. Lehnert, Strickland, & Armus 

2000). X-rav. JDahlem. Weaver. & HeckmaM Il998t 
Martin. Kobulnickv, & Heckman 2002 j), a n d ultraviolet 
wavel engths ([Heckman & LeithereJ 119971 iPettini. et. alJ 
2000). However, comparatively little is known about the 
eventual fate of the wind material — whether it escapes from 
the galaxy potential, or cools and eventually rains back down 
on the disk. A number of arguments have been made to 
this effect. IHeckman. Lehnert. Strickland. & Armusl (2000) 
found the terminal velocity of starburst-driven winds to be 
300-900 km s" 1 , independent of the host galaxy properties, 
and invoked simple models of the gravitational potential to 
suggest that galactic outflows escape from the potential wells 
of dwarf galaxies, but not from more massive hosts. However 
these simple scaling arguments neglect the multi-phase nature 
of galactic winds, and the complicated interaction of the wind 
with the interstellar medium of th e disk and the halo. Drawing 
on the numerical simulations of iMac Low & Ferrara (1999) 
which attempt to includ e this complicated gastrophysics, 
Ferrara & Tolstov (2000) conclude that galaxies with gas 
masses greater than ~ 10 9 M Q do not experience outflows. 
In contrast, our results imply that galaxies with masses as 
high as 10 10 M Q still eject some fraction of their metals into 
the intergalactic medium. This result is in line with the recen t 
observational and theoretical work of Strickland et al. ( 2003). 
These authors suggest that blowout is possible from galaxies 
with masses as high as 10 11 M based on a model where 
supernovae in disk OB associations work together to power 
the outflow. 

Of course it is quite likely that even in the absence of winds 
chemical evolution does not proceed in the simple manner we 
have assumed. It is worth considering if other factors could 
explain the trends we observe. For example, another means 
of lowering the effective yield is the inflow of metal-poor gas. 
While inflow may occur in some galaxies (the Milky Way, for 
instance) we do not regard it as a viable explanation for the 
correlation between effective yield and baryonic mass. Metal- 
poor inflow acts to lower the metallicity and to raise the gas 
mass fraction, so the net effect on the effective yield is modest. 
For example, consider a galaxy presently undergoing closed 
box chemical evolution with a gas fraction of 0.3. If such a 
galaxy experiences an inflow of metal-free gas with a mass 
equal to its present gas mass, the metallicity is reduced by 
50%, the gas mass fraction is increased by 54%, and the effec- 
tive yield is decreased by 22% (Alog^yy = 0.1 dex). Thus, 
even dramatic inflo w events have o nly a moderate effect on 
the effective yield. iGarnettl |2002) pointed out that the low 
effective yields of the smallest dwarfs would require them to 
accrete as much as 80-90% of their gas at late times without 
experiencing much star formation. We therefore consider it 
unlikely that metal-poor inflow is solely responsible for the 
low yields we observe. 

Another factor that can lower the effective yield is a break- 
down of the instantaneous recycling approximation - the as- 
sumption that gas is either immediately enriched and returned 
to the ISM or locked up forever in low mass stars. While this 
assumption obviously does not hold in detail, the main prac- 
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tical concern is at late times when low mass stars return their 
comparatively un-enriched gas to the interstellar medium. If 
this were a chemically significant effect then we would expect 
the ef fective yield to d e crease a s a function of mean galaxy 
age. fKauffma nn et alJ (l2003hl see Fig. 1) demonstrated a 
tight correlation between stellar mass and the D„(4000) spec- 
tral index, which is a good tracer of the mean stellar age of a 
population. The correlation is in the sense that more massive 
galaxies appear to be older. Hence we would predict that y e ft 
should decrease with galaxy mass, just the opposite of what is 
seen. While gas recycling may indeed influence the detailed 
shape of the y e ff—Mb aryon relation, it is almost certainly not 
the main driver of the observed trend. 

It is our view that the strong positive correlation between 
effective yield and baryonic mass is most naturally explained 
by the increasing potential barrier which the metal-laden wind 
must overcome to achieve 'blow out' . While the correlation 
is not particularly tight (±0.15 dex), its very existence nev- 
ertheless implies two very interesting things. First, nearly 
all low mass galaxies and many high mass galaxies have 
experienced some sort of blow-out event. Since the en- 
ergy requirements for blow-out are fairly exacting, it seems 
likely that these events are associated with starburst activ- 
ity. This chain of reasoning leads us to suggest that the star 
formation histories of most galax ies may be bur s ty rather 
than c ontin uous. (However, s ee Skillman (2001), van Zee 
(2001), and iHunter & Gallag herl i 19851) for some opposing 
viewpoints.) Unfortunately too little is known about the ex- 
tended gaseous halos of normal galaxies and the conditions 
requ ired for blow-out to m ake a more quantitative statement 
(see Strickland et al. 2003). The second significant implica- 
tion is that blow-out events have an important chemical im- 
pact even when integrated over a galaxy's history of star for- 
mation. This suggests that blow-out is either a very frequent 
occurrence or very catastrophic. In the coming years numer- 
ical and semi-analytical models should help to address these 
questions. In the meantime we emphasize our empirical find- 
ings: galactic winds are ubiquitous and extremely effective in 
removing metals from galaxies. 

The quest to understand how the chemical properties of 
galaxies couple to their star formation histories has been given 
added impetus of late by measurements of the metallicity- 
luminosity relation of galaxies at i ntermedi ate (0.3 < 
z < 1 .0) (Kobu lnickv et all 120031 Lilly. Carollo. & Stocktonl 
l2003tlMaier. Meisenheimer. & Hippelein|2004l) and high (z > 
2) redshifts llKobulnickv & Kool l2000t IPettini et all 120011) . 
Our results imply that metallicity is not a straightforward met- 
ric of galaxy evolution because metals can escape galactic po- 
tential wells. However, the strong correspondence of metal 
loss with the size of the potential has interesting implications 
if effective yields can be measured reliably. The combination 
of the intermediate and high redshift data with the correlations 
we measure at z ~ 0.1 will provide an important benchmark 
for successful models of feedback and galaxy evolution. 
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TABLE 1 

THE g-BAND Luminosity-Metallicity 
Relation 



M g 12 + log(Q/H) 





Pi <; 

1 Z.3 


Pis 

' ID 


1 Dli 


PHA 


Pai <; 


-22.37 


8.86 


8.99 


9.08 


9.15 


9.26 


-21.99 


8.88 


9.00 


9.09 


9.15 


9.24 


-21.61 


8.87 


8.99 


9.09 


9.15 


9.24 


-21.23 


8.84 


8.97 


9.07 


9.14 


9.23 


-20.84 


8.77 


8.92 


9.04 


9.13 


9.21 


-20.45 


8.71 


8.87 


9.01 


9.11 


9.19 


-20.05 


8.66 


8.82 


8.97 


9.09 


9.17 


-19.67 


8.60 


8.76 


8.93 


9.07 


9.17 


-19.26 


8.55 


8.71 


8.89 


9.05 


9.16 


-18.87 


8.39 


8.64 


8.83 


9.02 


9.14 


-18.46 


8.32 


8.61 


8.81 


9.00 


9.14 


-18.07 


8.21 


8.52 


8.70 


8.96 


9.11 


-17.68 


8.20 


8.41 


8.66 


8.89 


9.12 


-17.28 


8.18 


8.34 


8.62 


8.87 


9.07 


-16.87 


8.08 


8.23 


8.49 


8.81 


9.15 



NOTE. — The columns labeled P represent the 
2.5, 16, 50, 84, and 97.5 percentile of the distribu- 
tion of metallicity in bins of 0.4 dex in luminosity. 
Pso is the median. The g-band luminosity is in AB 
magnitudes and has been corrected for inclination- 
dependent extinction but not internal extinction. 



TABLE 2 

Thez-band Luminosity-Metallicity 
Relation 



M c ° r 12 + log(Q/H) 





Pis 


Pl6 


Pso 


Pu 


P97.5 


-23.67 


9.00 


9.06 


9.12 


9.18 


9.29 


-23.30 


8.95 


9.05 


9.11 


9.18 


9.27 


-22.91 


8.94 


9.03 


9.11 


9.17 


9.25 


-22.53 


8.92 


9.01 


9.09 


9.15 


9.24 


-22.13 


8.87 


8.98 


9.07 


9.14 


9.22 


-21.74 


8.83 


8.94 


9.04 


9.12 


9.19 


-21.35 


8.75 


8.88 


9.00 


9.09 


9.17 


-20.94 


8.70 


8.83 


8.96 


9.06 


9.15 


-20.56 


8.64 


8.77 


8.90 


9.02 


9.12 


-20.16 


8.58 


8.70 


8.85 


8.98 


9.10 


-19.76 


8.47 


8.64 


8.79 


8.94 


9.08 


-19.36 


8.37 


8.61 


8.73 


8.89 


9.06 


-18.96 


8.23 


8.52 


8.67 


8.83 


9.04 


-18.56 


8.15 


8.40 


8.62 


8.78 


9.02 


-18.18 


8.15 


8.31 


8.57 


8.72 


8.94 


-17.79 


8.16 


8.24 


8.45 


8.66 


8.87 


-17.39 


8.04 


8.15 


8.36 


8.59 


8.94 



Note. — The columns labeled P represent the 
2.5, 16, 50, 84, and 97.5 percentile of the distribu- 
tion of metallicity in bins of 0.4 dex in luminosity. 
P50 is the median. The ^-band luminosity is in AB 
magnitudes and has been corrected for internal ex- 
tinction. 



The Mass-Metallicity Relation 

TABLE 3 
The Mass-Metallicity Relation 



log(M»/M s ) 12 + log(Q/H) 





Pi 5 


p 16 


p 50 


p 84 


P ql 5 


8.57 


8.18 


8.25 


8.44 


8.64 


8.77 


8.67 


8.11 


8.28 


8.48 


8.65 


8.84 


8.76 


8.13 


8.32 


8.57 


8.70 


8.88 


8.86 


8.14 


8.37 


8.61 


8.73 


8.89 


8.96 


8.21 


8.46 


8.63 


8.75 


8.95 


9.06 


8.26 


8.56 


8.66 


8.82 


8.97 


9.16 


8.37 


8.59 


8.68 


8.82 


8.95 


9.26 


8.39 


8.60 


8.71 


8.86 


9.04 


9.36 


8.46 


8.63 


8.74 


8.88 


9.03 


9.46 


8.53 


8.66 


8.78 


8.92 


9.07 


9.57 


8.59 


8.69 


8.82 


8.94 


9.08 


9.66 


8.60 


8.72 


8.84 


8.96 


9.09 


9.76 


8.63 


8.76 


8.87 


8.99 


9.10 


9.86 


8.67 


8.80 


8.90 


9.01 


9.12 


9.96 


8.71 


8.83 


8.94 


9.05 


9.14 


10.06 


8.74 


8.85 


8.97 


9.06 


9.15 


10.16 


8.77 


8.88 


8.99 


9.09 


9.16 


10.26 


8.80 


8.92 


9.01 


9.10 


9.17 


10.36 


8.82 


8.94 


9.03 


9.11 


9.18 


10.46 


8.85 


8.96 


9.05 


9.12 


9.21 


10.56 


8.87 


8.98 


9.07 


9.14 


9.21 


10.66 


8.89 


9.00 


9.08 


9.15 


9.23 


10.76 


8.91 


9.01 


9.09 


9.15 


9.24 


10.86 


8.93 


9.02 


9.10 


9.16 


9.25 


10.95 


8.93 


9.03 


9.11 


9.17 


9.26 


11.05 


8.92 


9.03 


9.11 


9.17 


9.27 


11.15 


8.94 


9.04 


9.12 


9.18 


9.29 


11.25 


8.93 


9.03 


9.12 


9.18 


9.29 



NOTE. — The columns labeled P represent the 2.5, 16, 50, 
84, and 97.5 percentile of the distribution of metallicity in bins 
of 0. 1 dex in stellar mass. P 50 is the median. 
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TABLE 4 

The Baryonic Mass-Effective Yield Relation 



log <M baryBn M e ) logyeff 





' 2.5 


Pit. 


P™ 


Pa/i 


p 

'97.5 


8.67 


-2.84 


-2.74 


-2.57 


-2.34 


-2.08 


8.76 


-2.75 


-2.69 


-2.52 


-2.32 


-2.05 


8.86 


-2.80 


-2.68 


-2.51 


-2.27 


-2.03 


8.97 


-2.74 


-2.60 


-2.39 


-2.18 


-1.98 


9.06 


-2.72 


-2.60 


-2.40 


-2.21 


-1.98 


9.16 


-2.73 


-2.53 


-2.35 


-2.17 


-1.98 


9.26 


-2.69 


-2.49 


-2.31 


-2.13 


-1.92 


9.36 


-2.65 


-2.46 


-2.30 


-2.12 


-1.93 


9.47 


-2.61 


-2.43 


-2.27 


-2.11 


-1.95 


9.56 


-2.56 


-2.39 


-2.24 


-2.08 


-1.89 


9.66 


-2.55 


-2.38 


-2.24 


-2.09 


-1.90 


9.76 


-2.49 


-2.35 


-2.21 


-2.06 


-1.90 


9.86 


-2.50 


-2.33 


-2.19 


-2.04 


-1.87 


9.96 


-2.46 


-2.32 


-2.18 


-2.03 


-1.87 


10.06 


-2.44 


-2.30 


-2.17 


-2.02 


-1.85 


10.16 


-2.44 


-2.29 


-2.15 


-2.01 


-1.83 


10.26 


-2.43 


-2.27 


-2.13 


-1.99 


-1.82 


10.36 


-2.43 


-2.26 


-2.12 


-1.98 


-1.80 


10.46 


-2.41 


-2.25 


-2.11 


-1.96 


-1.80 


10.56 


-2.38 


-2.23 


-2.09 


-1.94 


-1.78 


10.66 


-2.38 


-2.22 


-2.07 


-1.93 


-1.77 


10.76 


-2.36 


-2.21 


-2.07 


-1.92 


-1.78 


10.86 


-2.36 


-2.21 


-2.06 


-1.91 


-1.77 


10.95 


-2.36 


-2.20 


-2.05 


-1.90 


-1.75 


11.05 


-2.38 


-2.20 


-2.05 


-1.90 


-1.74 


11.15 


-2.36 


-2.19 


-2.04 


-1.88 


-1.70 


11.25 


-2.38 


-2.20 


-2.03 


-1.86 


-1.73 


11.35 


-2.30 


-2.18 


-2.01 


-1.82 


-1.68 



NOTE. — The columns labeled P represent the 2.5, 16, 50, 84, and 
97.5 percentile of the distribution of the effective yield in bins of 0. 1 dex 
in baryonic mass. P 50 is the median. 



